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1 Introduction 

The one-dimensional Bose gas with delta function interaction (contact interac- 
tion) is a paradigmatic solvable [14] many-body quantum system. It shares its R- 
matrix underlying the integrability with the spin-| Heisenberg chain, but may be 
considered even simpler, since its Bethe ansatz equations have only real roots [10]. 
This feature, the absence of strings, in first place simplifies the thermodynamic 
Bethe ansatz (TBA) analysis, which is the reason why the Bose gas was the first 
non-trivial Bethe ansatz solvable model whose thermodynamics was analyzed in 
detail [18]. 

When in the 80s and early 90s the algebraic Bethe ansatz was applied to 
the study of correlation functions of solvable models, the Bose gas was again in 
the center of interest [10]. Explicit results, generalizing early work [12,13,15], 
were obtained in particular for the case of infinite repulsion, the so called impen- 
etrable Boson limit. In this limit many correlation functions can be expressed 
as Fredholm determinants or Fredholm minors related to integral operators of a 
special type, the so-called integrable integral operators [4,10,11], which have a 
close connection with classical integrable evolution equations and for this reason 
are suitable for an explicit calculation of the asymptotics of correlation functions. 

Another type of expression for the correlation functions of solvable many-body 
systems are the so-called multiple integral representations first obtained in [5] for 
the density matrix of a segment of the XXZ Heisenberg chain by an approach 
based on the representation theory of the quantum affine algebra Uq{sl2)- Later 
this result was rederived by means of the algebraic Bethe ansatz [9] which was 
important, since the algebraic Bethe ansatz turned out to be flexible enough 
for a number of generalizations. In our context we would like to mention the 
works [8] and [2], where multiple integral representations for a one-parameter 
generating function of the zz-correlation functions of the XXZ chain were derived 
for the ground state and for finite temperature. From [8] it was only a small 
step to the corresponding expression for the Bose gas [6], namely a multiple 
integral representation for a generating function of the density-density correlation 
functions. A finite temperature version was obtained directly from [2] in a certain 
scaling limit close to the ferromagnetic point A = — 1 of the XXZ chain [16]. 

The latter is the starting point of this work and is reviewed in the next section. 
Then, in section 121 we perform the impenetrable Boson limit and show that the 
resulting expression is equal to a known Fredholm determinant representation. 
This is the main result we wish to convey: in an appropriate limit the multiple 
integrals turn into a Fredholm determinant. Looking at it from a different angle 
this means, that we may interpret the multiple integrals as a 'deformation of a 
Fredholm determinant'. In section H] we complement our work with an account of 
the free Boson limit which is less trivial as might appear at first sight. In fact we 
derive a multiple integral formula with permanents replacing the determinants 
of the impenetrable case, which we were not able to spot in the literature. To 
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round off this work we demonstrate in section [5] how to obtain the density-density 
correlation function of the Bose gas from the generating function in the two 
hmiting cases of free and of impenetrable Bosons. Section [6] is devoted to a 
concluding summary. 



2 Density-Density Correlations in the Bose Gas 

The one- dimensional Bose gas with contact interaction is described by the Hamil- 
tonian 



H 



(1) 



Here ^/'^(z) and tp{z) are Bose fields with canonical equal-time commutation re- 
lations which act on the interval [0, i] for which we assume periodic boundary 
conditions, c > is the coupling constant. 
The operator 



Qix) 



dz ip^{z)ip{z) 



(2) 



measures the number of particles in the interval [0,x], < x < i. We need it 
to define a one-parameter generating function of the density-density correlation 
functions by {e'^^^^^)T,ij,, where the brackets indicate the grand-canonical ensem- 
ble average for a heat and particle bath of temperature T and with chemical 
potential /z. The function (e'^'^*^^^)^,/^ is particularly convenient in the context of 
the algebraic Bethe ansatz [10]. With the shorthand notation j{x) = ip'^{x)ip{x) 
for the particle density operator we have the following formulae, 

1 



i/j=0 



(j(O)j(^)) 



</3=0 



(3) 



In our previous work [16] we obtained the multiple integral representation 



oo ^ - n 

5^(dj2 n 

n=0 ^ ' 4=1 



27r 1 + e^fe)/^ 



27r 



R+iO 



n 



Vj - <ik 



ic 



det [M{pj, qk)] ^ det [G{pj, qu 



j,k=l,...n 



(4) 



for the generating function in the thermodynamic limit by considering a special 
scaling limit of the XXZ chain close to the isotropic ferromagnetic point A = — 1. 
Here M(p, q) is defined a^ 



M(p, q) 



(p-q) 



ic 



n 



p- qi + ic 



+ 



ice^ 



p - q + ic -^J^ p - pi + ic p-q 



n 



p-qi 



IC 



IC p-pi 



IC 



(5) 



"'^Note a typo in equation (26) of the published version of [16]: the first factor of the second 
term on the rhs should read ce"^ /{wj ~Pk){uJj —Pk ~ ic) instead of ce"^ /{wj ~pk){wj ^Pk +ic). 
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The temperature and the chemical potential enter through the functions e{p) 
and G{p, q) which must be calculated as solutions of integral equations. The 
'dressed energy function' e{p) is the solution of the non-linear integral equation 

dg c 



e[p) =p - - T 



ln(l + 



-e{g)/T 



) 



(6) 



of Yang and Yang [18]. The 'density function' G{p,q) solves the linear integral 
equation 

c f dk c G{k, q) 



G{p, q) 



{p - q){p - q - ic) 



+ 



TT 



(p - A;)2 + c2 1 + 



(7) 



3 Impenetrable Boson Limit 

It is known since the work of Girardeau [1] that in the limit c ^ oo the wave 
functions of Bosons with contact interaction turn into those of free Fermions, up 
to a function which takes values ±1. We therefore expect the multiple integral 
formula (jlj) to simplify in this limit. In the following we indicate the limit by 
supplying a superscript (oo) to the respective functions. 

We first of all note that the kernel in the integral equations ([6]) and ([7]) for 
e{j)) and G{j), q) vanishes for c — oo, such that their solutions become explicit, 



ei^){p)=p^-^^ G^^\p,q) 



—1 



p-q 



Consequently the expression 1/(1 + e^*^^-*/"^) turns into the Fermi function 

for non-relativistic particles. Moreover, the function M(p, q) simplifies drastically 



for c 



cxo, 



M^^\p,q) 



i(l -e^) 
p-q 



(10) 



and the exphcit product of the right hand side of equation (jlj) converges to one. 

Substituting ([8]) and (ITO!) into the determinants in (jlj) we see that they are of 
Cauchy-type in the limit. Using the well-known formula 

_ Ua<biPa-Pb){qb-qa) 



det 



Pj - Ik 



j,k=l,...,n 



Ulb=liPc^ - (lb) 



(11) 



we obtain 



IT, 11 

oo 

E 

ra=0 



;i -e^)" 
(n!)2 



n 

•i=i 



^e'^.-fip 



dq 
2^ 



A^(p) A2(g) 

Y{lk=M3-Pkf 



(12) 
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where we introduced the notation 



Hp) = lliPk 

j<k 



det[pi ],fc=l,...,n 



(13) 



for Vandermonde determinants. 

Due to the symmetry of the integrand with respect to all qj , one of the Van- 
dermonde determinants A(g) in each term can be replaced by a product of the 
diagonal elements of the Vandermonde matrix, and subsequently the q integrals 
can be pulled into the second Vandermonde determinant o Then 



1 



n 

■1=1 



dqi 
2tt 



-iqix 



A2(g) 



M+iO 



uik=iiQj-Pky 



det 



M+iO 



dq e-'^i^qi+k~^ 



(14) 



j,k=l,...,n 



Obviously the integrals inside the determinant on the right hand side can now 
be calculated by means of the residue theorem. Yet, it turns out to be more 
convenient to perform a number of elementary row- and column operations first, 
resulting in the sequence of identities 



det 



M+iO 



dq e-^^^qi+k-'i 

nr=i(?-p/)' 



j,k=l,...,n 



det 



dg e YliJi iQ-Pa) ULi il - Ph) 



27r 



A2(p) 
1 

A2(p) 



+iO 

det 



]Xi=i{<i-piy 



j,k=l,...,n 



dq 



-iqx 



R+iO 



2vr {q - Pj){q - Pk) 



det 



2 sin 



XI 



j,k=l,...,n 



Pj - Pk 



(15) 



j,k=l,...,n 



In the last line for j = A; it is understood to take the analytic continuation of the 
function 



V{u,v) 



(16) 



U — V 

namely V{u, u) = x, which comes from the calculation of the residua at the 
second order poles of the diagonal entries. Finally, pulling out the exponential 
factors and a minus sign in f[T^ . we end up with the expression 



det 



R+iO 



dg e-'''^g^+'=-2 

ULiiq-pi)' 



j,k=l,...,n 



-l)«A2(p) 



det [V{pj,Pk)].,^^^ 



(17) 



^The same trick was used for the free Fermion limit of the XXZ chain in [7]. 
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We further introduce the notation 



Vf{u,v) = y^)V{u,v)^W) (18) 
and substitute it together with and (fTTI) into (fT2|) to obtain 



lT,lJL 



n=0 



n 

.j=i 



27r 



det [VF(Pj,Pfc)],.fc^i 



(19) 



Now recall (see e.g. [17]) that for an integral operator K with kernel K{p,q) 
acting on a function defined on an interval / as 



im{p)= jdqK{p,q)<P{q) 



(20) 



its Fredholm determinant has the infinite series representation 

det(l + ir) = V- / dpi...dp„det[ir(p„pfc)]^^^^ ^. (21) 
n=o ^' 

Then 



(e-'3(-))g = det(l + ^V» 



(22) 



where Vp is the integral operator with kernel Vf(p, g). This is a known result [10]. 
What we find remarkable, however, is that it directly follows from the multiple 
integral representation (jl]). 



4 Free Boson Limit 

We first of all note that the limit c — in (jlj) exists. In the following we indicate 
it by a superscript (0) to the respective functions. Observing that, in this limit, 
the kernel in the integral equations ([6]) and ([7j) for e{p) and G{p, q) turns into a 
representation of the 5-function, their solutions become explicit. 



r£(o)(p)i 






[ T \ 


= exp 


[ T \ 



1 + exp 

and the matrix simplifies, 

M(o)(p„gfc) 



G(o)(p,g) 



— c 



[p - g)2 1 - e-(p'-^'yT 



(23) 



1 -e^ 
(Pj - Ik)' 



Ua=liPj-Qa) 



Pa) 



(24) 
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Substituting (123!) and into the integral representation (jl]) for c ^ and 
using (fTTj) we obtain 



E 

n=0 



1 -e^)" 
(n!)2 



dp. 



n ' 

det 



27r e(P?-/^)/^ - 1 
1 



27r 



'IQjX 



det 



jr',A:=l,...,n 



1 -2 



(25) 



j,k=l,...,n 



To simplify the integrand under the q integrals we utilize an identity borrowed 
from \8], 



'- a<b 



'-a,b=l 



Y[ (Pa -Pb + ic) sign((T) J(p^(j) , Qk) 



Yl{Pa-qb){Pa-qb + ic) 



det 



a,b=l 



.{Pj - Qk)iPj - Qk + ic) 



, (26) 



j,k=l,...,n 



where the summation extends over all permutations a G &n and where 

IILl [ UtZliPa -qb + ic) Ub=a+liPa " Qb) 



HPj^Qk) 



Ua<biPb - Pa + ic) 



(27) 



The arguments pj, qu are understood to represent the sets {pj}-, {qu}- Performing 
the free Boson limit c — > in fl2B]) we obtain a useful expression for the quotient 
of determinants under the integral in 



det 



1 



.{Pj-flkf 



det 



j,fc=l,...,n 



Pj - qk 



E 



1 



j,k=l,...,n CTge^ 



riLib-w -?a) ■ 



(28) 



Squaring of (l28i) produces two independent sums over permutations. Then, 
proceeding similarly as in the case c — oo, we may replace 



E 

crgSn 



E 



]Ta=l{Pa - Qa) ]Ia=l{P.(a) ' Qa) 



(29) 



under the multiple-g integrals due to their symmetry. Consequentially, the q in- 
tegrals factorize and can be evaluated by means of the residue theorem as in 

(USD, 



dg 



-iqx 



27r {u — q){v ~ q) 



-r(M,T;)e-'("+'')"/^ 



(30) 
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Finally, introducing the shorthand notation 

1 



b{p) 



(31) 



the multiple integral representation of the one-parameter generating function of 
the correlation function in the free Boson limit reads 



yQ{x)\(0) 



E 

n=0 



1)" 



.i=i 



perm[VB(pj,pfc 



j,k=l,...,n ■ 



(32) 



Here the chemical potential must be restricted to the physical range /i < for the 
integrals involving the Bose function b{p) to converge. Equation (1321) looks rather 
similar to our former expression f|T9|) for impenetrable Bosons. The determinants 
are replaced by permanent^ and the Fermi function is replaced by the Bose 
function here. 



5 One- and Two-Point Functions 

In this section we would like to recall briefly how to obtain the one- and two-point 
density correlation functions by applying ([3]) to f|T9l) and fl32|) . Because ip appears 
therein only in the factor (e*^ — 1)", the derivatives with respect to if can be easily 
calculated by means of the formula 



a^(e^-i)" 



6n,i, dl{e^-ir =Sn,l + 26n,2. (33) 

ip=0 



Then we obtain for the impenetrable Boson case the known [10] explicit expression 
(J(^))f:^ = 5./ ^fiP> = I ^f{p) = D,{T,f,), (34) 

which due to translational invariance is independent of the interval length x. The 
analogous well-known expression for free Bosons, 

(J(^))S = / ^b{p) = Ds{T,f^), (35) 
can be obtained from (l32l) by application of ([3]). 



•^The permanent of an n x n matrix A = (^"' j,)j,fc=i,...,n is defined by 
E ^'-(i) ■ • -^"-(n) = perm [^'' ■ 
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As the term arising from n = 1 is linear in x we need to consider only the 
n = 2 term in order to calculate the density- density correlation function. For the 
impenetrable Boson case we have 



dpi 


f dP2 


R 






f dpi 




' 271 



fiPl) f{P2) 



2 



{.P1-P2Y _ 



dp 
2^ 



fip) 



2 / ^ /" ^ f{pi) fip2) COs{piX - P2X) 



2D|(T,/i)-2 



dp 
2^ 



f{p) cos{px) 



(36) 



implying another well-known result [10], 



0-(o)j(^))S-^^(T,/i) 



dp 
2^ 



f{p) e 



ipx 



(37) 



for the density- density correlations. Similarly, the corresponding expression for 
free Bosons can be obtained from ([3]) and fl32|) . 



(j(o)j(^))?i-^l(r,/i) = + 



dp 
2^ 



b{p) e 



ipx 



(38) 



6 Conclusion 

In the limit c ^ 00 of impenetrable Bosons we transformed a multiple integral 
representation for a certain generating function of the density-density correlation 
functions of the Bose gas with contact interaction into a Fredholm determinant of 
an integrable integral operator. The latter provides a link to classical integrable 
evolution equations and thus a means for a rigorous analysis of the asymptotics of 
correlation functions [3]. Taking this into account, it is tempting to interpret the 
multiple integral representation as a 'deformation of a Fredholm determinant'. 
One may hope that such an interpretation will finally help to solve the problem 
of the calculation of the long-distance asymptotics also in the generic case of 
finite repulsion. 

For the free Bose gas, c 0, we rearranged the summands in the multiple 
integral representation in terms of permanents. Interestingly the resulting series 
has some similarity with the Fredholm determinant of the impenetrable case. 
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The Fermi function is replaced by the Bose function, and the determinants are 
replaced by permanents. 

In the generic case of finite c > our original integral representation (jlj) 
interpolates between the two limiting cases f|T9|) and fl32|) . What is special about 
the limiting cases is, that only the first two terms in the infinite series contribute 
to the density-density correlation function. This property seems lost in the generic 
case. We expect, however, that taking into account only a few terms of the series 
may provide good approximations for small and for strong finite coupling. 
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